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We report our study on unconventional charge density waves (UCDW) (i.e. a charge density 
wave with wavevector dependent gap) in pure quasi-one dimensional conductors. We develop a 
new possible mechanism of establishment of such a low temperature phase, in which the driving 
force of the phase transition is the electron-phonon interaction with coupling depending on both 
the momentum transfer (q) and the momentum of the scattered electron (k). Mean field treatment 
is applied to obtain the excitation spectrum, correlation functions such as the density correlator 
and the optical conductivity, and the effective mass of the phase excitation. The fluctuation of the 
order parameter leads to the sliding of the UCDW as a whole. In the absence of impurities, we 
calculated the effect of this fluctuation on the optical properties. The inclusion of the collective mode 
significantly alters the optical conductivity, and leads to an effective mass which is nonmonotonic 
in temperature as opposed to conventional CDWs. 

PACS numbers: 71.45.Lr, 73.20.Mf, 78.30.-j 

I. INTRODUCTION 

As a result of intense research during the past few decades, much is known about the properties of density 
waves (DW) possessing a constant A order parameter^ On the other hand, it is also well known that 
the momentum dependent order parameter A(k) plays an important role in theories of superconductivity 
(termed unconventional)^ Indeed, many unusual features of the high temperature superconductors are 
understood in terms of a k-dependent single particle gap of d-wave symmetry, and therefore the case of d-wave 
superconductors for high-T c cuprates is now well established 4£ Also, most heavy fermion superconductors 
and organic superconductors appear to be unconventional or nodal with gapless quasiparticle spectrum, 6 

It follows from the success of these developments naturally, that investigating the density wave sector 
may as well turn out to be a fruitful enterprise. Momentum dependent order parameter in an electron-hole 
condensate was first introduced in the context of excitonic insulator^ Since then considerable attention has 
been focused on investigating unconventional density waves (UDW) under various circumstancesASiifiiiAiiS 
A striking feature of these unconventional condensates is that due to the vanishing momentum average 
of the gap over the Fermi surface ((A(k)) = 0), these systems are not characterised by spatially periodic 
modulation of either the charge or the spin density, though the name "density wave" is widely used^ This 
property makes UDW a very likely candidate for those systems where a clear and robust thermodynamic 
phase transition is not accompanied by an order parameter detectable by conventional means, the situation 
often referred to as " hidden order" H 

Possible materials with UDW ground state include: the quasi-two dimensional transition metal dichalco- 
genid 2H-TaSe2 for which an /-wave UCDW was proposed 1 ^; the mysterious low temperature micromagnetic 
phase of the heavy fermion compound URu2Si2 for which an USDW has been suggested originally 15,16 then 
later further corroborated^; and the pseudogap phases of (TaSe^Ir^*^ an d the underdoped high-T c su- 
perconductors for which the d-w&ve DW scenario was proposed^. 

Perhaps, one of the most likely candidate for possessing some kind of CDW of unconventional type in the 
low temperature phase is the organic conductor cv-(BEDT-TTF)2KHg(SCN)4. This salt has been investigated 
experimentally extensively in the past few yearsrSiSi* 2 ^^ 2 ^!^ 5 . and recently, the experimental findings of these 
latter works have been explained rather well by a quasi-one dimensional UCDW with momentum dependent 
order parameter, including the threshold electric field, 26 ' 27,28 the angular dependent magnetoresistance ) 2 ^' 36 
and the magnetothermopower and Nernst effect^i. 

In this paper, our goal is to develop the theory of a novel type of quasi-one dimensional UCDW, in which 
the phase transition is governed solely by electron-phonon interaction. A preliminary report of some of our 
results have already been presented^ The article is organised as follows: in section 2 we define our model 
and develop its mean field theory and thermodynamics. Section 3 is devoted to the optical conductivity and 
the effective mass of the phase excitation. Our conclusions are given in section 4. 
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II. THE MODEL 

A. The coupling and the order parameter 

Let us start by considering the Hamiltonian of the interacting electron-phonon system of our quasi-one 
dimensional conductor 

i? = 5Z £ ( k ) c t<T C k^ + H LJA ( c l) a q,A a q.A + -^ ^(k,q)c+ +qCT c k , ff (a q . A + a+ q J, (1) 

k.cr q,A k,q,A,er 

where the one particle energy e(k) for the electrons is given by the usual highly anisotropic (t a ^> tf, ^> t c ) 
nearest neighbor tight-binding formula measured from the chemical potential 

e(k) = — 2t a cos(ak x ) — 2% cos(bk y ) — 2t c cos(ck z ) — [i. (2) 

Henceforth, the most conducting crystal axis x will often be called the chain direction. Our system is based 
on an orthorhombic lattice with lattice constants a, b, c, and with one atom per unit cell. The three acoustic 
phonon branches that the crystal possesses, labelled by the polarization index A, have all been taken into 
account in Eq. Q). As usual, Ck,er is the annihilation operator of a Bloch electron with momentum k and 
spin a in a single band, while a qjA is the same for a free phonon with energy uj\(q), momentum q and 
polarization A. The interaction matrix element in Eq. Q is given by 

-j=DA(k,q) = 1 / d 3 r^ +q (r)e A (q)hVL/ at (r)]Vk(r), (3) 

A A /2Mw A (q) J 

with [7 a t being the atomic potential, N is the number of unit cells in the sample, and e A (q) denotes the 
unit vector of polarization A and momentum q. The k-dependence of the coupling - emanating from second 
quantization on Bloch-basis (^k(r)) - turns out to be crucial in order to form an unconventional charge 
density wave (UCDW). Indeed, it is clear from Eq. J2J), that choosing plane waves (~ e lkr ) for the one 
electron basis, one ends up with the well known Frochlich model 33 - 34 of the conventional charge density 
waves (CDW) with coupling simplified to the Fourier component (at wavenumber q) of the atomic potential. 
This approximation would therefore lack any kind of unconventionality, which can only originate from the 
k-dependence of the electron-phonon coupling. 

Let us now proceed with the mean field treatment of the Hamiltonian in Eq. . We introduce the phononic 
operator (f>\ (q) = + a+ q A . Within mean field theory, the operator is replaced by its expectation value 
(</> A (q)), which differs from zero only if the momentum is set to q = ±Q, with Q = (2kp, ir/b, tt/c) being 
the best nesting vector. With this, the electronic part of Eq. simplifies to a quadratic Hamiltonian in 
the electronic operators, and one obtains 

ffMF = ]T *+(k)(C(k)p 3 + A'(k) Pl - A"(k)p 2 )f (k). (4) 
k 

Here pi stand for the Pauli matrices acting on the space of left-, and right-going electrons, \I/ + (k) = 
( c k-p C k-Qt> c ki' C k-Qj.) ^ s ^ ne ^ our component spinor operator, £(k) denotes the linearized electron spec- 
trum in k x around ±fcp, the Fermi surface: £(k) = ^fCMeI — kp) — 2tbCOs(bk y ) — 2i c cos(cfc z ), and finally 
A(k) = A'(k) + iA"(k) is the decomposition of the order parameter into real and imaginary parts. The 
order parameter itself reads as 

A(k) = J= 5^ A (Q))^Ck, -Q). (5) 

Now the diagonalization of Eq. (@J is straigthforward: for the excitation energies we have £7±(k) = ±£'(k), 
E(k) = -\/£(k) 2 + |A(k)| 2 , which is the usual two band quasiparticle spectrum known for example from 
the theory of unconventional spin-, and charge density waves driven by purely electronic correlations.^ 
The thermodynamics of the system in the low temperature phase is determined solely by the momentum 
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dependence of E(k), and the possible node structure (line or point nodes) on the Fermi surface, if there are 
any. Thus in the followings we explore the momentum dependence of the order parameter. Making use of 
the Bloch- Wannier transformation for the electronic wave functions, ipk(r) — 1/VN J2 r e lkR </?(r — R), we 
change representation in Eq. J3J to Wannier basis and obtain 

^(k,-Q)=g A 5> A (R)e-* kR , (6) 



R 



where 



A A (R) = e A (Q) J d 3 r^(r) (^T e lQR ' [-Vf/ at (r - R')]j p*(r-R), (7) 

and g\ = (2Mw A (Q)) _1 / 2 . Equation © is clearly a Fourier expansion of the coupling in its argument k 
in the Brillouin zone harmonics, with coefficients being overlap integrals of Wannier orbitals localized on 
different sites. Since our system is a tight-binding solid (see Eq. J2Jl), these are considered to be exponentially 
small if the spacing is large between the two orbitals. Therefore it is a reasonable approximation to retain 
only the on-site (R = 0) and the nearest neighbor terms (R = ±a, ±b, ±c). Indeed, as is easily seen from 
Eqs. (SJ and © , in order to keep the possibility of an UCDW formation with momentum dependent order 
parameter, one has to go beyond the on-site term. This situation is similar to the development of a CDW in 
an interacting electron system: in the simple Hubbard model the interaction U < is not able to drive the 
system into an unconventional charge density wave ground state either, and one has to include two center 
exchange integrals as well to facilitate an unconventional condensate Aiii^ With all this 

D* x (k, -Q) = D$ + D$ cos(bk y ) + D\ sm(bk y ) + D$ cos(ck z ) + D% sin(cfc z ), (8) 

where we have set k x — kp, since the relevant k-dcpcndcncc of both the coupling and the gap is confined to 
a narrow region near the Fermi sheet at +/sf- The coefficients are therefore given by 

= gx[A x (Q) + A x ( a )e- lk - a + A x (-a)e^ a ], (9) 

D$ = g x [A x (b) + A x (-b)}, (10) 

D x 2 = ~igx{A x (b) - A x (-b)}, (11) 

D$ = g x [Ax(c) + A x (-c)], (12) 

D% = -ig x [Ax(c) - A x (-c)}. (13) 

Consequently, the gap will be of the form 

A(k) = A + Ai cos(bk y ) + A 2 sin(Wc y ) + A 3 cos(cfc z ) + A 4 sin(cfc z ), (14) 



A ' = ^E^(Q))A A , i = 0,...,4. (15) 



with 



N 

The structure of the gap function in Eq. I|14|) makes it clear, that there are line nodes on the Fermi surface 
in any of the latter four unconventional cases (e.g.: A2sin(6fc y ) = 0, k z varies freely). The quasiparticles 
around these nodes, often termed nodal in the terminology of superconductivity, will determine the nature of 
the thermodynamics. Since the gap function itself is the same as the one in quasi-one dimensional electronic 
UDWs- (either spin or charge), the thermodynamics of a one-component (with only one type of gap amplitude 
being finite among the four possible unconventional cases in Eq. (|14|l . e.g.: A(k) = A2sin(6fcj,)) phononic 
UCDW is therefore identical to those of either an electronic UDW, 9 or a d-wave superconductor— (in spite 
of the different topology of their Fermi surfaces). 
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B. The gap equation 

In this subsection, we shall now proceed with the derivation of the self-consistency condition for the order 
parameter A(k) known as the gap equation. Let us consider the equation of motion of the phonon operator 
0a(Q) introduced in the previous subsection 

-^^a(Q)=^a(Q) 2 0a(Q) + 2^ a (Q) j B a (-Q), A = 1,2, 3 (16) 

with 

Bx(-Q) = -!= J2 ^(k, -Q)4_ Q , CT c k , CT . (17) 

k,(T 

Taking the expectation value of Eq. 1(16(1 over the full Hilbert space, the anomalous electron-hole expectation 
value appearing in Eq. 1(17(1 can be easily calculated in the low temperature phase, and using Eq. J^J) we get 



(18) 



At this point it is important to call the attention to the fact that since the susceptibility matrix xw m 
Eq. (|18fl is clearly not diagonal, the three equations for the corresponding phonon modes are thus coupled 
together in Eq. 1)16(1 . This can be easily understood: although there is no direct interaction between the 
different phonon modes, the electron system due to its coupling to all branches acts as a medium, and 
mediates an effective phonon-phonon interaction. Applying the ansatz (4>\{t)) = (4>\(0))e~ int to Eq. (|16|) 
averaged over the full Hilbert space, the eigenfrequencies (n = 1, 2, 3) of the new, effectively noninteracting 
quasiparticles are obtained from the following eigenvalue problem 

" 2 W>a) =XK*aa'+2u; a Xaa']<<M, A = 1,2, 3. (19) 

A' 

Due to the temperature dependence of xaa' , the matrix in the square brackets in Eq. 1(19(1 may develop a zero 
eigenvalue signalling a soft mode at sufficiently low T. The point where this takes place defines the transition 
temperature T c of the charge density wave instability. The gap equation is essentially equivalent to the 
softening condition of any of the three new phonon modes below the transition temperature (£!(T < T c ) = 0) 

= ]TK<5aa'+2xaa'](<M, A = 1,2, 3. (20) 

A' 

Adding up the three equations for A in Eq. I(2U(I . after some straightforward algebra one obtains the gap 
equation for A(k) in the familiar forrrASi 

A(k) = lVP(k,kOA(kV anh( ^ (k/)/2) , (21) 
v ) N v i ; v ; 2£(k / ) ' v ' 

P(k, k') = Y, HW^(k, -Q)D x (k\ -Q). (22) 
Once the gap equation is solved, comparing Eqs. 1(21122(1 with Eq. JSJ yields the soft mode amplitude as 
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At this point, it is worth noting, that one can deduce the very same gap equation by minimizing the 
thermodynamic potential of the system 



with respect to its undetermined variables (<fi\). 

Before we go on to the next subsection dealing with the solution of the gap equation, we find it necessary 
and instructive to summarize what the theory has provided us with so far. Up to this point our calculations 
have been completely general in the sense, that no assumptions were made on the particular form of the 
electron-phonon coupling D\ (k, q) . The primary formula in Eq. J5j) shows that the momentum dependence 
of the gap stems solely from the coupling. If it is unconventional, then the averaged value over the Fermi 
surface vanishes ((A(k)) = 0) leading to a spatially homogeneous density of the electronic charge. This state 
of affairs might remind one to the similar result found in unconventional spin-, and charge density waves 
governed purely by electronic correlations^ Although there are similarities in many aspects between the 
two types of condensates, there is one point that is completely different in nature in the phononic UCDWs 
under focus, and must be emphasized by all means. Namely, though there is no spatial modulation of the 
electronic density at all, the Peierls distortion of the underlying lattice with periodicity Q is still present 
(note the finite soft mode amplitude in the low temperature phase in Eq. 123(1 ). which could be detected by 
usual means, like x-ray scattering. This feature is to be contrasted with that of an electronic DW (either 
spin or charge), where the density wave instability is caused purely by the electron-electron interaction, and 
the lattice is not affected at all. Further, if the system favours an unconventional ground state, then there 
is no modulation either in spin or charge. Thus these systems are often referred to as "systems with hidden 
order" , and are consequently promising candidates for systems where robust thermodynamic signals of a 
phase transition are seen without any order parameter. 

On the other hand, the static deformation of the lattice is a charasterictic feature of the conventional 
charge density waves with constant gapp and in such systems it is always accompanied by the electronic 
charge oscillation. Now we face a different and novel situation, namely: in a phononic UCDW these two 
phenomena need not go together. 

Regarding all these, we conclude that the notion "hidden order" might not be applicable to the phononic 
UCDW, because as entering the low temperature phase a new periodicity with wavevector Q develops. 
Nevertheless, the model exhibits all peculiar properties of the aforementioned UCDW that were applied to 
explain the experimental findings measured on the a-(BEDT-TTF) 2 KHg(SCN) 4 salti 26 i 27 i 28 i 29 i 30 i 31 More- 
over, the extra feature of the distorted lattice might give an explanation to a recent x-ray study performed 
on this salt, 25 in which enhanced structural modulation is found below T ~ I0K suggesting a coupling to 
the electronic degrees of freedom. 

In the followings, the next subsection will be devoted to the solution of Eq. J2U- 



In this subsection we focus on the solutions of the gap equation obtained in Eq. 121fl . The knowledge of the 
explicit momentum and temperature dependence of A is essential for our further investigations. So far the 
calculation is completely general within the frames of the Hamiltonian in Eq. For further study however, 
we retain only the longitudinal term in the interaction in Eq. 1|22[1 . which is considered to be the strongest 
component, and ignore the coupling to the transverse phonon modes. Therefore we drop the polarization 
index A in the followings for brevity. With this, the kernel simplifies to P(k, k') = 4lu(Q)~ 1 D* (k)D(k'). 
The immediate consequence of the factorization of P in its momentum arguments, that follows right from 
the gap equation and from the explicit form of D*(k) (see Eq. (JSJ) is 



^,T;{(^)}) = i^c A (Q)|<^)| 2 -| Yl Ml + e- a/3£(k) ) 




(24) 



C. Competing phases in the ground state 



Aj(r) = a 



i,j = 0,...,A. 



(25) 



The relation in Eq. 1251) not only fixes the ratios of the different gap amplitudes to a given, temperature 
independent value, but also involves that there is only one T c , where each type of phases (with nonzero 
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FIG. 1: Left panel: the temperature dependence of the function h(T; {Do, D2}) for fixed DijTc — 0.7 and for 
different D0/D2 values, where ho(T) = h(T; {0, -D2}). For the numerical computations hf^f/Tc = 30 was applied. 
Right panel: the quasiparticle density of states for A0/A2 = 0, 0.6, and 1.3. The curve with the highest peak at 
\u>\ = A2 belongs to the value (one component UCDW), while the one with symmetrically placed smaller peaks 
belongs to 0.6, finally the curve with one peak and a clean gap belongs to 1.3. 



coupling) simultaneously open. Accordingly, the solution of Eq. l(2"B|) is simply Ai(T) = Dih. Here h(T; {Di}) 
is a dimensionless function of order unity determined by 



from which the transition temperature is obtained as 

Tc = ?l VF k F e-" (Q)/PM0)<|B| 2 )]. (27) 

In equation H27JI Po(0) = a/irvF stands for the normal state density of states per spin, and (...) denotes the 
Fermi surface average. Furthermore, at zero temperature and around T c we have 

7r / (\D\*hx(T-i\D\)) \ 
h(0) = - exp ( 71Tn ^ ) , (28) 




hlT) " T imim^-T- <29 » 

The whole temperature dependence of h below T c is shown in Fig. ^ (left panel) for the case where only 
Do and D2 are set to nonzero, real and positive values in Eq. (JHJ. This restricted choice of the parameters 
is sufficient and simple enough to present results for a coexisting CDW+UCDW phase. Consequently, the 
order parameter will have the form 

A(k) = A + A 2 aiii(6Ay) 1 (30) 

with the gap amplitudes both positive. At this point it is important to call the attention to the fact, that 
since the coupling amplitudes Di in Eq. JSJ) can be complex as well, it is even possible that there is a relative 
phase between A and A 2 leading to an optical gap in the single particle density of states regardless of the 
values of |A | and |A 2 |, but this case is beyond our scope of investigation in this paper. We note however, 
that complexity of the same order parameter in electronic UCDW can account for a number of peculiar 
properties observed in (TaSe^I^ 
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Since A(k) = h(T)D*(k), thus Eq. © implies that (4>) = \fNh(T) for the longitudinal phonons. Conse- 
quently the Q Fourier component of the lattice distortion can be written as u(Q) ~ g((f>) e(Q) ~ h(T). This 
Peierls distortion of the ions can be seen by x-rays even if the electron system is unconventional. 

Now we calculate the quasiparticle density of states for the coexisting phase having the type of gap function 
in Eq. (|3L)fl . The density of states reads as 



Po(0) Jo 2^ |A + A 2 sm(y)\< 




(31) 



where K(z) is the complete elliptic integral of the first kind. The energy dependence of p{uj) is shown in 
Fig. ^ (right panel) for different values of Ao / A 2 . It is clear from the figure and can be easily verified from 
Eq. (|31|1 . that a true optical gap G — 2(Ao — A 2 ) opens at the Fermi energy only if Ao > A 2 , furthermore 
the positions of the logarithmically divergent peaks are given by 



| Wpcak | = 9(A 2 - A )|A ± A 2 | + 6(A - A 2 )(A + A 2 ). (32) 

One can also easily check, that in the Ao = limit (Do = 0), Eq. (|31|l simplifies to the well known result 
of a single component unconventional density wave^ and obviously the A 2 = limit (_D 2 = 0) yields the 
conventional BCS result. 



III. CORRELATION FUNCTIONS 



In this section we calculate the frequency dependent conductivity of the previously introduced model, 
namely that of the coexisting CDW+UCDW phase with A(k) = Ao + A 2 sin(6fcj,). The outline is what 
follows: we shall start with the formulation of the Green's function. We then proceed with the evaluation of 
the quasiparticle contribution of the density-density correlator. We do this, because we are primarily focused 
on the optical conductivity in the chain direction, which can be determined from the charge conservation 
rule, that is the continuity equation. Afterwards we study the one bubble result for the conductivity, we go 
on and incorporate the effect of interaction between electrons and phonons into the theory on the level of 
random phase approximation (RPA). The RPA calculation enables us to identify the effective mass of the 
collective motion, the sliding of the density wave. 



A. Single particle conductivity 



The Green's function of a charge density wave, either conventional or any type of unconventional, is given 
by G _1 (k, iujn) — iu n — H(k), where H(k) = ^(k) j 03 + A'(k)/9i — A"(k) j o 2 is the four by four matrix appearing 
in the mean field Hamiltonian in Eq. Q). Here, unconventionality is manifested solely in the momentum 
dependence of the gap. From this, G is explicitly obtained as 

rf , ■ , + g(k) P3 + A'(k)pi - A"(k) P2 

^n + £(k) + |A(k)| 2 

Now, the single particle contribution to the retarded product of the density correlator, IFq = ([n, n])o, is 
calculated in the usual way and reads as^ 

11°! (C, W -» W ) = -j^ Tr ( G ( k ' lLU n)G(k + q, lLU n + IV)) 

= 2p (0) 7 2-^(l-/), 

C — UJ Z 
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where, for simplicity, we limited our analysis to q = (q x ,0,0) (i.e. wave vector pointing in the quasi-one 
dimensional direction) and C — vfQx- In addition, / is the generalized version of the function that also 
appears in the correlation functions of conventional DWs with constant gap^S only minor differences are 
present due to incorporation of unconventionality 



/ = (C 2 -^)-/ / tanh -— -Re ' yy ' = dydE, (35) 



Wo io V 2 J D % fW- 

N = (C 2 -to 2 ) 2 -4E 2 (( 2 +w 2 ) + 4C 2 |A(y)| 2 , 
D = N 2 — 64E 2 lu 2 ( 2 (E 2 — |A(y)| 2 ), 



where A(y) = A + A 2 sin(?/). 

Now, making use of Eq. 1)34(1 and charge conservation, the single particle contribution to the complex 
conductivity in the chain direction is of the form a(^,oj) = ne 2 m~ 1 ioj(u! 2 — £ ) (1 — /), where n/m — 
2po(0)vp/V c , with n being the particle density and V c is the cell volume. Furthermore, in the £ — ► long 
wavelength limit, the real and imaginary parts of the optical conductivity (cr = a% + io-i) are given by 

01 = D5(w) + {ne 2 /m)^- 1 ^, (36) 
a 2 = (ne 2 /m)uj- 1 (l-.Q, (37) 

where /o = lim£_>Q / — /o + */o' ^ s the decomposition of /o into real and imaginary parts, and D = 
(ne 2 ir/m)[l — /o(0)] is the Drude weight. The regular part of <j\ in Eq. (|36l) can be calculated explicitly, 
since /q has a closed form if both gap amplitudes are real. That is, we again constrain our analysis to a real 
gap function, just as we did in the previous section (see Eq. I|30|l and subsequent argument). With all this 

/o'H = tanh(^)M^)i. (38 ) 

Here 



C(x,u) = —= = \2[{u + iy + x 2 ]K 2,. 

^/(x+i) 2 -u 2 1 1 v y (x+i) 2 -u 



+ I) 2 - AE (2^/=^=) - Mu - X + 1)U (-1^,2^=^=:) } , (39) 

x = |o;|/(2A2), u — A0/A2, K(z) and E{z) are the complete elliptic integrals of the first and second kind, 
while II(n, z) stands for the incomplete elliptic integral of the third kind. The imaginary part ai in Eq. I|37|l is 
evaluated numerically with Kramers-Kronig transformation at zero temperature. The frequency dependence 
of both quantities are shown in Fig. [3 for different values of the ratio A0/A2, including the case of Ao = 0, 
which is simply the one component UCDW limit It is interesting to see from Fig. [21 that whatever small 
the A2 is, while we are in the Ao > A2 regime, the sharp onset of absorption with square root divergence at 
the gap edge known for a conventional density wave is disappeared. The size of the gap is G — 2(Ao — A2), 
that has been mentioned already along with the single particle density of states (see Fig.^). In addition, the 
peaks of logarithmic type are apparently located at 2w poa k, where w pea k is the corresponding singular point 
in the quasiparticle density of states defined in Eq. (|32|l . 



B. Conductivity with collective contribution 



We shall turn now our attention to the effect of collective contributions on the results of the previous 
subsection. Namely, we investigate how the conductivity is affected if we include the short wavelength 
component of the interaction within RPA. This part of the interaction is responsible for the density wave 
instability and depending on the strength of the couplings in Eq. (JSJ), it drives the system into one of the 
possible symmetry breaking ground states, either a conventional CDW, a single component unconventional 
CDW, or even one of the many coexisting phases. We proceed with the assumption we made through the 
paper, that is the ground state of the interacting electron-phonon system is a coexisting CDW+UCDW 
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FIG. 2: Left panel: the real part of the optical conductivity in the quasi-one dimensional direction at zero temperature 
for A0/A2 = 0, 0.6 and 1.3, respectively. The curve with a single peak at \u\ = 2A2 belongs to the value 0, and is 
the result for a single component UCDW. The curve with two symmetrically placed peaks around 2A2 belongs to 
0.6. Finally the spectrum with a clean gap around zero frequency has Ao = I.3A2. Right panel: the imaginary part 
of the optical conductivity at zero temperature for the same Ao / A 2 values that appear on the left panel. The line 
with a positive jump of tv/2 at \u\ — 2A2 is the single component UCDW result with Ao = 0. The other two belong 
to 0.6 and 1.3 (from left to right). 



phase with order parameter A(k) = Ao + A2 sm(bk y ). Our goal is to calculate the dressed chain direction 
conductivity by taking into account the fluctuation of the order parameter, leading to the sliding of the 
condensate as a whole. With this, we are in a position to identify the effective mass of the collective 
excitation, the same way as it was done in the seminal paper of Lee, Rice and Anderson. 34 

The interaction part of the Hamiltonian in Eq. Q , describing scatterings of the electrons from one Fermi 
sheet to the other, can be recast as 



H ' = -7^ E^ k ' -Q)* + ( k + q)p-*(k)0+(Q) + h.c.}, 

q,k 



(40) 



where p± = (pi ± ip^jl, and only the coupling to the longitudinal phonon has been retained in accordance 
with what has been pointed out in the subsection dealing with the solutions of the gap equation. Namely, 
it is a reasonable approximation to omit the coupling to the transverse modes, as they do not really form 
ionic charge fluctuations that the electrons could feel and couple to. The bare phonon propagator, that acts 
as the interaction between electrons in the diagrammatic language, is given by 

X>°(q ) tw n ) = - / dT e l ^(T T 0(q,T)0(-q)) 



2w(q) 







With all these, the coupled RPA equations for the renormalized density-density correlator II11 
are 



(41) 
([n,n]) 

(42) 
(43) 
(44) 

Here LT^ is given by Eq. I|34l) . and the other zeroth order correlators (with superscript 0) are also evaluated 
within the standard method, and we find 



n u = - n° £>°n_i - n°_2?°n +1 , 



nn 


= n 


n+i 


= n 


n_i 


= n 



n° 



P U IE 



n_i = n°_i - n°_ + x>°n_i - n°_£> n +1 . 



n°t T (c,w 



1 



|£(k)| 2 Tr(p ± G(k,z W „)p T G(k + q,z, 



iv)) 



^(Q) Po(0) 
2 4/i 2 



(45) 



[(- 2 -C 2 )/"2/ 4 ] 
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M/A 2 




|w|/A 2 



(d) 
















|w|/A 2 |w|/A 2 

FIG. 3: The real (left panels) and imaginary parts (right panels) of the optical conductivity in the RPA at zero 
temperature from Eqs. 15011 and I5H , The A0/A2 ratio is fixed to the same values as in Fig. [5] (one bubble results), 
that are 0, 0.6, and 1.3. Beside that, the dimensionless, renormalized coupling appearing in the RPA formulae is set 
to A(0) = 0.4 for figures (a) and (b), A(0) = 0.7 for (c) and (d), and A(0) = 1.1 for (e) and (f), respectively. 



and similarly 



n; ±(C >^) = ±^c/. 



(46) 
(47) 
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1.2 









D /D 2 = 0.6 
















J 0.2 


0.4 0.6 0.8 




T/T c 



FIG. 4: . The effective mass versus reduced temperature for both a conventional CDW and a CDW+UCDW phase 
with A(0) = 0.4. In the latter case, we present results for A (T)/A 2 (T) = D /D 2 = 0, 0.6, and 1.3 with D 2 /T c = 0.7 
fixed. The CDW result is plotted for comparison, where we used D 2 = 0, and Do/T c = 0.7. For each curve 
m*(0)/ro = 3.5. 



In addition, f± differs from / only in the integrand in Eq. (|35[) . namely there is an extra |A(y)| 2 factor in 
the numerator. After some algebra, one finds for the dressed density correlator 

n " = 2 ^c^ i 1 - ^ira ) ' (48) 

where A = Ao[w 2 (Q)/(|2A| 2 )] = w(Q)/9o(0)/(2/i 2 ) is the renormalized temperature dependent electron- 
phonon coupling, while A = 2p (0)(\D\ 2 ) /lu(Q) is the bare coupling, see Eq. (|2Tjl . With the aid of Eq. (gHJ) 
one ends up with the following formula for the inchain optical conductivity 

a{o J ) = —-(-l + h-h^), (49) 
where m* /m = 1 + A -1 /^ 1 is the effective mass. The real and imaginary parts (cr = ci + ia%) read as 

" = ™<"> + I ^HTTw- (50 » 
, 2 = ^I( 1 _4±W), (51) 

m uj \ |1 + A/or / 

with £) = (ne 2 7r/m)[l — X~ 1 m/m*] the Drudc weight. Both <j\ and are shown in Fig. Inversus frequency 
at zero temperature, for different values of A(0) and A0/A2. Now we have to point out one interesting 
feature of <Ti, that cannot be seen properly from the figure. As the function /q has logarithmic singularities 
at frequencies u> = 2oj pca k, furthermore /q is bounded, this implies that the regular part of <j\ in Eq. (|5UI) 
vanishes in an inverse logarithmic manner at the same frequencies. Since the logarithmic singularity is very 
weak, the plots in Fig. cannot really resolve this kind of behavior around the zeros. Nevertheless, as the 
coupling A(0) increases, the feature becomes more and more apparent. We shall note at this point, that 
the way how the RPA result relates to the quasiparticle contribution shown in Fig. |2 can be contrasted 
to what has been found for a conventional charge density wave^i In such a system, the familiar square 
root singularity at the gap edge 2A is suppressed and transformed to a square root edge by the collective 
contributions. In our phononic UCDW however, the absorption is finite for the whole energy range, unless 
Ao > A2, where a clean optical gap develops in the spectrum. 

The effective mass of the sliding for our CDW+UCDW model can be seen in Fig. 0] with A(0) = 0.4. The 
parameters used for the numerical calculations are: Aq(T)/A2(T) = D0/D2 — 0, 0.6, and 1.3, respectively 
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with D%/T c = 0.7 being fixed (see the results for h(T/T c ) on Fig.^ left panel). Along with these new results, 
we have also plotted the effective mass of a conventional CDW for comparison, where we used Dq/T c = 0.7. 
Looking at the figure, an interesting new feature of the phononic UCDW can be established: as long as Dq 
is smaller than D2, that is the system possesses small energy excitations, Dirac fermions, around the nodes 
of A(k), the effective mass m* is a nonmonotonic function of the temperature. That is not the case for a 
conventional CDW, as it decreases in a monotonic fashion as we approach T c from below4£ 

IV. CONCLUSIONS 

In this paper we have developed the mean field theory of a novel type of quasi-one dimensional uncon- 
ventional charge density wave, where the density wave instability is caused purely by phonons. As the 
electron-phonon coupling exhibits significant dependence on the momentum of the scattered electron, this 
leads us to a wavevector dependent order parameter taking different values on different points on the Fermi 
surface. In particular, it turns out that the explicit form of the single particle gap is identical to that in an 
electronic UCDW, where in contrary, the phonons do not play any role in the development of the density 
wave condensate, leaving the lattice unaffected. It is important to realize, that if the system favours an 
unconventional ordering in the ground state, then due to the vanishing momentum average of the order 
parameter the system lacks any periodic oscillation in the electronic charge. Though the two type of con- 
densates have apparently much in common, like thermodynamics determined by the Dirac electrons around 
the nodes of the gap for instance, there is an important difference that has to be emphasized. That is in a 
phononic UCDW, while the electronic density remains homogeneous, the underlying ionic lattice undergoes 
a distortion as we enter the low temperature phase, clearly signalling the phase transition. For this reason 
we might say that the "hidden order" reveals itself in such a way, and could be experimentally accessible 
and measurable by x-ray scattering. 

On the other hand, the aforementioned Peierls distortion is a key feature of the conventional charge density 
waves arising from the coupling of the lattice to the electronic degrees of freedom. Hence, we conclude that 
the present model of the phononic UCDW is a natural generalization of the conventional theory to more 
complicated, even k-dependent couplings. 

We have calculated the optical conductivity in the quasi-one dimensional direction, and identified the 
effective mass of the collective phase excitation. The calculations were carried out using a multicomponent 
gap allowing us to investigate the effect of the coexisting conventional and unconventional orders. Con- 
sequently, the spectra are more structured compared to the pioneer result of Lee, Rice and Anderson^ 
Nevertheless, though the singularities and zeros are of logarithmic type characteristic to UDWs in general, 
the single component UCDW limit shows the familiar property: the singularity at the maximum optical gap 
is suppressed by the collective contributions. In addition there is considerable absorption for small frequen- 
cies arising from the nodal excitations. We have also found that the effective mass exhibits a novel type of 
nonmonotonic temperature dependence as long as the unconventional order dominates the low temperature 
phase, as opposed to conventional systems. This could serve as a valuable tool in identifying the nature of 
the low temperature phase. 
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